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Abstract. 

The existence of a scaling evolution for cosmic string loops in an expanding universe 
is demonstrated for the first time by means of numerical simulations. In contrast 
with what is usually assumed, this result does not rely on any gravitational back 
reaction effect and has been observed for loops as small as a few thousandths the size 
of the horizon. We give the energy and number densities of expected cosmic string 
loops in both the radiation and matter eras. Moreover, we quantify previous claims 
on the influence of the network initial conditions and the formation of numerically 
unresolved loops by showing that they only concern a transient relaxation regime. 
Some cosmological consequences are discussed. 



PACS numbers: 04.50.+h, ll.lO.Kk, 98.80.Cq 
Cosmic strings are line-like topological defects which may have formed during 
phase transitions in the early universe. These defects are primordial vacuum remnants 
associated with the spontaneous breakdown of particle physics symmetries occurring 
during the cooling of the universe [H [2|. If cosmic strings are formed at the end 
of inflation [3], they should still be present in the observable universe and many 
studies have been devoted to the understanding of their cosmological implications (see 
Refs. in |5] for a review and references therein and thereto). Recently, cosmic strings 
have been of particular interest since it was pointed out that Fundamental and Dirichlet 
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strings formed at the end of brane inflation could also play a cosmic string-like role in 
cosmology O [71 [8] . 

It is now part of the common lore that the existence of a string network is 
cosmologically acceptable owing to the so-called "scaling regime" of "long strings". 
Intersections between super-horizon sized strings (the long, or infinite, strings) produce 
sub-horizon sized loops, ensuring that the total energy density of long strings scales with 
the cosmic time t as instead of the naively expected and catastrophic a being 
the scale factor. Nevertheless, overdosing of the universe is avoided only if the energy 
density in the form of loops is radiated away. This is typically the case for oscillating 
loops by means of gravitational wave emission but may no longer be true if the strings 
carry particle currents [9l [TOl [TT|, [12] . 

On sub-horizon sizes, the understanding of string networks is still a matter of 
debate. Early analytical studies predicted the scaling property of long strings by 
considering the network dominated by one length scale only: the inter-string distance, 
somehow related to the horizon size [T3l [H] . However, the first high resolution numerical 
studies identified a copious production of small sized loops, suggesting that the network 
dynamics should also involve a small fundamental length scale [I5l [161 [13 UHl [19] . 
This led to the development of the so-called three scales model [20]. In addition 
to two horizon sized distances, the long strings develop, through their intersections, 
a small scale wiggliness structure, explaining the formation of the small sized loops. 
An interesting feature of the model is that the small length scale reaches a scaling 
regime only if gravitational radiation is considered, otherwise, and under reasonable 
assumptions, the kinky structure is expected to incessantly grow with respect to the 
horizon size. The main features of the three-scales model have been numerically 
confirmed in Minkowski space-time [21]. However, in these simulations, nearly all loops 
are produced at the lattice spacing size, which makes the evolution and the scaling 
properties of the small scale structure strongly dependent on the cutoff. It was thus 
pointed out that if this feature persists whatever the lattice spacing, the typical size of 
physical loops might be the string width, which implies that particle production rather 
than gravitational radiation would be the dominant mode of energy dissipation from a 
string network [2T1[22]. Although loop formation at the string width has been recovered 
in more recent field numerical simulations, it has been argued that loop production 
rather than particle emission should remain the dominant cosmological energy carrying 
mechanism [23l [21] . 

In the following, we present new results coming from numerical simulations of 
cosmic strings evolution in a Friedman-Lemaitre-Robertson-Walker (FLRW) universe. 
For this purpose, we used an improved version of the D. Bennett and F. Bouchet 
Nambu-Goto cosmic string code [18]. We show for the first time evidence of a scaling 
regime for the cosmic string loops in both the radiation and matter eras down to a 
few thousandths of the horizon size. The loops scaling evolution is similar to the long 
strings one and does not rely on any gravitational back reaction effect. It only appears 
after a relaxation period which is driven by a transient overproduction of loops, with 
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Figure 1. The (lOO^c)^ comoving volume in the matter era when the observable 
universe occupies one eighth of the box. 

respect to the scaling value, whose length is close to the initial correlation length of the 
string network. We then discuss the effects induced by the finite numerical resolution 
and show that they do not affect the loops scaling regime. Moreover we confirm an 
explosive-like formation of very small sized and numerically unresolved loops during the 
first stage of the simulations, suggesting that particle production may briefly dominate 
the physical evolution of a string network soon after its formation. 

Our numerical simulations of strings in FLRW space-time are performed in a fixed 
unity comoving volume with periodic boundary conditions. The initial scale factor is 
normalised to unity while the initial horizon size is a free parameter which controls the 
starting string energy within a horizon volume. During the computations, the comoving 
horizon size grows and the evolution is stopped before it fills the whole unit volume for 
which the finiteness of the numerical box starts to be felt. We used the Vachaspati- 
Vilenkin (VV) initial conditions where the long strings path is essentially a random 
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walk of correlation length £c, together with a random transverse velocity component of 

root mean squared amplitude 0.1 |25l |18]. Such a velocity component allows a faster 

relaxation of the string network toward its stable cosmological configuration. The results 

presented below mostly come from two high resolution runs in the matter and radiation 

eras, respectively, performed in a (100£c)^ comoving box and with an initial string 

sampling of 20 points per correlation length (ppcl). The initial size of the horizon is 

chosen to reduce the relaxation time toward the scaling regime and has been set equal 

to c/ho = 0.063 for the matter era and (iho = 0.057 for the radiation era. The dynamic 

range achieved by these runs is 8 and 17 in conformal time, which is equivalent to 520 

and 308 in physical time, for the matter and radiation era, respectively (see Fig. [1]). 

In Fig. [2], we have plotted the evolution of the energy density p^o associated with 

the long strings, defined to be super-horizon sized, together with the loop energy density 

distribution dpo/da. More precisely, dpo(a) is the energy density carried by the loops 

whose length, expressed in units of the horizon size, (ih, is in the range a to a + da. 

A logarithmic binning in a of resolution Aa/a ~ 10~^ has been used in the range 

[10^^, 10^] to compute these quantities. From Fig. [2l the scaling values for the energy 

density associated with the infinite strings are 

£■ (fi 
p^^ = 28.4 ±0.9, Poo 77 = 37.8 ±1.7, (1) 

U mat (J rad 

for the matter and radiation era runs, where U denotes the string mass per unit 
length. The central value quoted above has been defined to be the mean value in 
the last half conformal dynamic range of each run while the errors correspond to the 
dispersions around the mean during that time. These scaling values are compatible 
with previous results while being slightly lower in the radiation era, possibly due to our 
better numerical accuracy and dynamic range [HI [12]. The really new results concern 
the loop energy density distribution which, after a transient regime (the bump in Fig. [2]), 
reaches a self-similar evolution where (\po{c^)d\ remains stationary, for all values of a 
down to approximately 5 x 10~^. The existence of such a loops scaling regime implies 
that 

^ = 5(a)-^, ^ = ^, (2) 
da ^ ' dl da arf^ ' ^ > 

where dn(a)/da is the number density distribution of loops. For a flat FLRW space- 
time, recall that cih = 3t in the matter era and = 2t in the radiation era. Notice 
that the loops smaller than 5 x 10~^ the size of the horizon have not yet entered such a 
scaling behaviour at the end of the runs. As can be seen in Fig. [21 the transient regime 
characterised by the bump in the energy density evolution is longer for the smaller loops. 
However, it is reasonable to believe that these loops would have reached the scaling 
regime with more dynamical range. These relaxation effects and their relationship with 
the initial conditions are discussed later in the letter. 

In order to determine the scaling function S{a) we have plotted in the top part 
of Fig. [3] the rescaled distribution arfj^dn/da as a function of a, and at equally spaced 
physical times spread over the dynamic range of the simulations. As expected from the 
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Figure 2. Evolution in the matter and radiation eras of the energy density associated 
with long strings and the energy density of loops of physical size /phys = Q^h- The 
time variable is the rescaled conformal time r?/^c- The energy density of long strings 
rapidly reaches a scaling regime where poa l/dl^, up to some damped relaxation 
oscillations. These plots show that the energy density of Zphys-sized loops also reaches 
a long string-like scaling regime where dpo oc This loop scaling regime appears 

after a relaxation period during which the energy density increases and decreases (the 
bumps) . This transient energy excess signs the relaxation of the initial string network 
towards its cosmological stable configuration. Note that the transient regime is longer 
for the smaller loops, as expected from the hierarchical process of loop formation. 
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establishment of scaling observed in Fig. [21 the distribution functions in Fig. [3] start to 
superimpose at the largest length scales during earlier times. Simultaneously, the non- 
scaling parts of the distribution function shift toward small a. This essentially means 
that the intersection processes which create loops from long strings take place in a very 
similar manner for loops themselves. Their self-intersections give rise to more numerous 
smaller loops in such a way that a constant energy flow cascades from the long strings 
to the smallest loops. Conversely, due to their small size, the reconnection processes of 
these freshly formed smaller loops to long strings and large loops is negligible [H] . As a 
result, since the source of energy is the entering of long strings inside the horizon size, it 
is not really surprising that the horizon ends up being the relevant quantity also for the 
loops. Note that some hints of these effects have recently been observed on the kinky 
structure of strings in Minkowski space-time p6] . 

We have performed a power law least squares fit of ad^dri/da in the domain 
Us < a < ttoo where the loops scaling regime has been reached. Our prescription is to 
set fts as the lowest size of loops, in units of the horizon size, for which the corresponding 
energy density dpo(«s) remains stationary during the last 5% of the simulation conformal 
time range (see Fig. [2] and Fig. [3]). In order to avoid contamination by long string 
effects, the upper bound has been fixed to the typical distance between infinite strings, 
i.e. = (U/ pooY^'^ /dh. By fitting the data points averaged over the above-mentioned 
5% conformal time interval, one gets the power law fit S{a) = Co with 



for the matter and radiation run, respectively. The errors quoted above mainly come 
from systematic errors associated with the determination of the a domains where the 
loops are in the scaling regime. They have been estimated as the largest Co and p 
deviations between four power law fits. These other fits use the same data but involve 
the four a intervals which are obtained by multiplying and dividing each of the original 
bounds Os and Ooo by a factor of two. 

As can be seen in the top panels of Fig. [3l the transient overproduction of loops 
preceding the scaling (the bump right to maximum) and the overall maximum of the 
loop distribution evolve in time. During the runs, they peak at decreasing sizes with 
respect to the horizon size. In the bottom panels of Fig. [3], we have plotted the rescaled 
distributions C = a^^^df^dn/da for the matter and radiation eras. For the length scales 
in scaling, these distributions are approximated by the above-mentioned constant Co 
whereas the relaxation regime and other non-scaling features appear as deviations from 
this fiat distribution. Note that, this time, we have plotted the rescaled distributions 
with respect to the physical length /phys = ad^- They now peak around a constant value 
which is close to the initial physical correlation length = 1/100 associated with the 
VV initial conditions. The setting up of the scaling regime appears in this frame as 
the progressive damping of the bump, which flattens the distributions toward Co and 
smooths out the correlations coming from the initial conditions. Moreover, the overall 
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Figure 3. The rescaled loop number density distributions with respect to a (top 
panels) and Zphys = ckc'h (bottom panels). They are plotted for different times starting 
at i = 1.1 and t — 0.8 (short dashed curves), with a physical time sampling equal to 1.1 
and 0.8 for the matter and radiation runs, respectively. The scaling regime propagates 
from the large scales toward the small scales, while the relaxation bump around the 
initial correlation length is progressively damped. The a intervals used to determine 
the scaling function in Eq. ^ are represented, as well as the best power law fit (long 
dashed line) and its estimated systematic errors (dotted lines). Note that the overall 
maxima of the loop distributions in the top panels correspond to the knee centered 
around the initial resolution length in the bottom panels. 



maximum of the arfj^dn/da distribution appears as a knee, which also remains at a 
constant physical length during the simulation. This length ends up being close to the 
initial resolution length associated with the VV initial conditions (see also Fig. H]). As 
detailed in Ref. [18], our numerical simulations are not lattice size limited in physical 
space due to an adaptive griding on each loop. However, at any time, loops with 
strictly less than 3 points cannot be formed and this is especially relevant for the initial 
conditions. Indeed, for an initial sampling of 20 ppcl, this cutoff necessarily leads to 
additional correlations in the VV network around the physical length = (3/20)£c, 
which precisely corresponds to the knee of the rescaled loop distributions C. The fact 
that the correlations associated with both and 4 remain at constant physical lengths 




Figure 4. Influence of the initial resolution length £,■ on the rescaled loop distributions 
at the end of three small (AOic)^ radiation era runs having an initial sampling of 10, 
20 and 40 ppcl, and a dynamic range of 45 in physical time (top panels). As in the 
previous figures, a — Zphys/c^h is the loop length in units of the horizon size. In the 
bottom panels, the influence of the 3 points cutoff is quantified for the 20 ppcl run by 
checking for stress energy conservation. 



during the subsequent evolution may be compared to the evolution of the small-scale 
correlation length obtained in the three-scale models in absence of gravitational damping 
and under other reasonable assumptions [20] . 

Since the above-mentioned discretisation effects concern the smallest loops, they 
should not influence the properties of the string network on larger scales. In the top 
panels of Fig. HI we have plotted the rescaled loop distributions obtained at the end 
of three (40£c)^ radiation era runs having an initial sampling of 10, 20 and 40 ppcl, 
respectively. As expected, one can check that the finite resolution effects remain confined 
to length scales smaller than the initial correlation length ic of the string network and 
do not affect the previously discussed loop scaling regime. We have also checked for 
these smaller runs the insensitivity of the loop distribution with respect to the initial 
random velocity component. 

In previous studies of Nambu-Goto strings, the dissipation effects have been 
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associated with the formation of lattice sized loops in physical space pT| [22] . In the 
context of our 3 points cutoff, all triangle shaped loops are removed from the subsequent 
evolution since they cannot self-intersect and do not significantly alter the evolution of 
larger scales, but would otherwise require additional computing time. Since such a 
removal is not equivalent to a fixed physical size cutoff, another way to discuss the 
dissipation effects is to test the total stress energy conservation during the evolution. 
Designating by p and P the total string energy density and pressure, we have plotted 
in the bottom panels of Fig. H] the evolution of the total network mass AA. = a^p, 
the total pressure work W = a^P, and the energy dissipation rate Q = M.' + 3?iW 
for the (40£c)^ radiation era run (20 ppcl). The prime stands for the derivative with 
respect to the conformal time and Ti = a'/ a is the conformal Hubble parameter. From 
the conservation of the Nambu-Goto stress tensor, one would expect Q = 0. As a 
result, the existence of a sharp negative peak for Q at the very beginning of the run 
signals a strong energy loss rate in the form of numerically unresolved loops (see Fig. H]). 
However, this happens only during a brief period in which the universe expands less 
than a factor of 10^'^ and corresponds to the times at which the overall maximum 
of the loop distribution around the initial resolution length appears. As the pressure 
evolution shows, this observation confirms the existence of an explosive-like relaxation 
of the VV strings network toward the cosmological network. Since in the early universe 
a string network soon after its formation is certainly far from its cosmological stable 
configuration, one might also expect such an explosive relaxation to physically take 
place. Note that during the later times of the numerical simulation the stress energy 
conservation is quite well satisfied, especially during the loop scaling regime. 

In summary, we have numerically shown that the energy and number densities of 
cosmic string loops in an expanding universe reach a scaling regime down to a few 
thousandths of the horizon size. This result does not rely on gravitational back reaction 
effects and suggests that the loop formation mechanism is the dominant energy carrying 
mechanism in the cosmological context. The loops scaling regime takes place after a 
relaxation period during which we observe an overproduction of loops peaked around a 
constant physical length close to the initial correlation length of the string network. 
We have confirmed the existence of a brief explosive-like dissipative behaviour in the 
very first stage of the evolution where a few percent of the total energy density are lost 
in the form of numerically unresolved loops. At the same time, the overall maximum of 
the loop distribution appears around a constant physical length scale close to the initial 
resolution length of the string network. These effects suggest that particle production 
may dominate the relaxation of a string network in the early universe from the string 
forming phase transition toward its cosmological scaling configuration. We have however 
shown that these discretisation effects do not affect the loops scaling regime. 

Our results should provide an improved basis for the development of loop evolution 
models [27] and for the determination of the gravitational wave background associated 
with the existence of a cosmological string network [28] . On the other hand, our Nambu- 
Goto string model certainly breaks down for length scales close to string width. One 
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therefore expect the cosmological loop distribution given in Eq. ([2]) to be valid above a 
given cutoff fixed by physical processes not included in our analysis. 
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Note added 

A week after our results were rendered public on the archives, two other papers dealing 
with the scaling of cosmic string loops appeared: Refs [29] and [30] . Each of these works 
exhibits its own loop distribution and we would like to briefly discuss qualitatively what 
we suspect to be at the origin of these differences. 

In Ref. [29], the authors studied cosmic string evolution under the simplifying 
hypothesis of Minkowski spacetime and observed that most of their primary loops are 
produced at a size 0.3t before going into loops of scaling size O.lt, a value close to their 
long inter-string distance. In order to attempt a direct comparison between Minkowski 
and FLRW simulations, one may express lengths with respect to the typical inter-string 
distance (our ttoo): their primary loop production size 0.3t would correspond to our 
a ~ 0.5 (for the matter era). As can be seen in our Fig. [3l this is out of our scaling 
range. Nevertheless, making abstraction of the statistical fluctuations induced by the 
small number of loops we have on that scale, a peak might be guessed, at least for the 
matter era run. As far as we can compare these results, they seem to be associated 
with loop formation events typical of the one-scale Kibble model. On the other hand, 
their loop distribution at small scales differs from ours and we notably do not observe 
the equivalent of their peak around O.lt. In our opinion, this may be related to their 
numerical cutoff which switches off the loop dynamics in the network as soon as their 
size is smaller than 0.25t. We do not use any similar assumptioiii|. Let us also notice 
that their long dynamical range is obtained at the expense of a procedure to increase 
the simulation volume, a procedure which may introduce spurious correlations on all 
length scales and which is hardly tractable in FLRW spacetime. 

X More recently, Ref. fSTj presented some new results by using the FLRW integration scheme of Ref. [18] . 
Their loop distributions are still different from ours albeit closer. The Q.2ht cutoff has been replaced 
by a removal of all loops which do not intersect after one oscillation. 
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In Ref. |3n], the authors performed the numerical simulations in FLRW spacetime 
and their results are more directly comparable to ours. They used a high number of 
points per correlation length (ppcl) to sample the strings and therefore have a smaller 
dynamical range than us. As a result, the loops we have found to be in scaling might 
still be in the transient regime at the end of their runs (see our Fig. [2]). The number of 
ppcl has two effects. Firstly, high resolution ensures high numerical accuracy to describe 
the Nambu-Goto dynamics, i.e. the numerically computed worldsheet evolution remains 
closer to the true solution for a given algorithm. As thoroughly discussed in Ref. [18] 
(Sect. IIB-3), the integration method we used reaches the same numerical accuracy as 
the other codes but with a factor of three to four less points per correlation length. The 
second effect of large ppcl is to reduce the initial resolution length and consequently 
the typical size of the purely numerical extra-correlations added in the initial string 
network by the discretisation. There is no mention in Ref. [30] of such kind of effect. 
However, their Fig. 4 is almost identical to our loop distribution of Fig. [3] but plotted 
against l/t. From this figure, they claim that the peak location of their loop distribution 
is moving slower at late time as if the whole distribution were approaching a scaling 
regime. In our interpretation, since we have found that the overall peak is precisely 
located at 4, a fixed physical length, such a slowing down solely comes from the choice 
of variable C/t (see our Fig. [3l bottom panels). 

In spite of these differences, both of these works claimed to observe some scaling 
evolution for the cosmic string loops, independently of radiative back reaction effects, 
thereby confirming our main result. 
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